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摘 要 本 文通 过 SST 水流 模型 求解 RANS 方程 组 , 研究 了 带 径 向 轮 缘 密 封 的 预 旋 盘 腔 流 动 换 热 特性 , 并 采用 附加 变量 
法 研究 径 向 密封 的 封 严 特性 . 比较 了 不 同 冷 气量 , 不同 进 气 预 旋 比 下 的 盘 腔 内 流动 系数 、 转 盘 表面 努 塞 尔 数 和 封 严 效率 . 
结果 表明 : 在 中 分 面 处 ,， 预 旋 比 随 进 气 预 旋 比 的 增 大 而 增 大 ， 静 压 系数 沿 径 向 线性 增 大 . 主流 入 侵 使 转盘 壁面 绝热 温度 升 
高 , 根据 流动 状态 计算 得 出 的 Nw 较 实际 略 大 . 进 气 预 旋 提 高 封 严 效 率 . 
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Abstract 
numerically investigated by RANS equations and SST turbulence model, and the additional passive 


The flow and heat transfer in a pre-swirl rotor-stator cavity with rim radial seal are 
tracer method is used to simulate the mainstream ingestion phenomenon. The variations of flow 
coefficient in cavity, the Nusselt number on rotor surface and sealing efficiency with different sealing 
flow rate or pre-swirl ratio are investigated. The numerical results show that in the split, pre-swirl 
ratio increases with the increase of the ratio of inlet pre-swirl. Static pressure coefficient along the 
radial direction increases linearly. On account of the mainstream ingestion which makes adiabatic 
temperature on rotor rise, The nusselt number calculated by the flow state is slightly larger than the 
actual. Inlet pre-swirl is beneficial to sealing efficiency. 
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涡轮 盘 腔 冷 却 是 燃气 透 平 冷却 系统 的 重要 组 成 
部 分 . 在 实际 流动 中 ， 主 流 的 高 温 燃 气 经 由 轮 缘 处 
的 密封 结构 入 侵 轮 盘 腔 室 ， 同时 主流 燃气 对 轮 盘 边 
缘 的 加 热 作用 会 使 得 轮 盘 温度 升 高 , 热 应 力 增 大 , 最 
终 导 致 轮 盘 变形 , 存在 安全 隐患. 在 盘 腔 冷却 中 , 最 
常用 的 方法 是 向 盘 腔 中 通 入 从 压气 机 中 引出 的 冷却 
空气 ,冷却 空气 进入 盘 腔 对 转盘 壁面 进行 冷却 ， 同 
时 作为 封 严 气体 对 轮 盘 腔 室 密 封 。 通过 减少 从 压气 
机 引入 的 冷气 量 ， 可 以 提高 透 平 气动 效率 ， 因 此 对 
盘 腔 内 流动 与 换 热 特性 的 研究 具有 重要 意义 。 

国内 外 学 者 已 经 对 盘 腔 的 流动 换 热 和 轮 缘 密 封 
的 封 严 效率 两 个 方面 进行 了 多 次 实验 及 数值 研究 ， 
在 轮 缘 密封 的 封 严 效率 方面 ，Patinios 等 中 采用 
低 雷 诺 数 相似 的 实验 装置 对 1.5 级 双重 密封 的 盘 腔 
系统 进行 了 封 严 效率 的 实验 研究 , 发 现 静 边界 层 内 气 
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流 轴 向 穿 过 核心 区 到 达 动 边界 层 ， 同 时 核心 区 预 旋 
比 决定 盘 腔 径 向 压力 分 布 ， 因 此 对 核心 区 的 流动 状 
态 研 究 十 分 重要 .Sangan 等 2-4 进行 了 单 级 透 平 
实验 ， 对 轮 缘 密封 外 部 诱导 入 侵 (ET)、 旋 转 诱 导入 
侵 (RD 和 双重 密封 的 封 严 特性 进行 了 实验 和 理论 分 
析 , 提出 了 经 实验 验证 的 孔 板 模型 . Paul 等 回 为 得 
到 高 质量 的 非 定常 轮 缘 密封 模型 ， 实 验 测 量 盘 腔 内 
外 压 差 随 密封 内 流速 的 定常 变化 ， 并 对 转动 时 盘 腔 
内 压力 分 布 频率 ,进行 频谱 分 析 。Svilen 等 加 通过 
实验 测量 , 比较 简单 密封 和 双重 密封 对 封 严 气量 、 叶 
片 压力 场 、 大 尺度 下 的 不 均匀 入 侵 、 非 设计 工 况 下 的 
主流 流动 变化 时 的 敏感 度 。 

在 流动 与 换 热 方面 , Liu 等 [1 对 预 旋 喷 嘴 位 于 盘 
腔 不 同 径 向 位 置 的 实际 绝热 系统 进行 研究 ， 其 预 旋 
温 降 与 功 耗 呈 线性 负 相 关 . Zhang 等 &9] 通过 CFD 
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数值 计算 ， 比 较 了 不 同 长 径 比 接收 孔 时 排 气 系数 、 
绝热 效率 等 流动 和 换 热 特性 的 变化 ， 研 究 了 实际 结 
构 的 盘 腔 系统 不 同 预 旋 角 度 下 盘 腔 的 整体 流动 情况 ， 
确定 预 旋 角度 为 20° 时 性 能 最 为 优越 . Tao 等 td 运 
用 PIV 实验 方法 和 CFD 数值 方法 , 研究 了 自由 转 
静 盘 腔 系 统 的 风阻 温 升 ， 研 究 得 出 转盘 温 升 受 盘 腔 
间隙 影响 较 小 ， 仅 为 在 一 定 清流 系数 下 ， 旋 转 雷 详 
数 和 通 流 雷 诺 数 的 函数 . Luo 等 9 采用 TLC 技术 
对 实际 结构 的 低位 进 气 预 旋 盘 腔 进行 实验 ， 研 究 了 
不 同 转速 与 冷气 量 对 盘 腔 换 热 特性 的 影响 ， 给 出 了 
实验 得 到 的 Nwu 转盘 分 布 云图 .Owen 等 03 改进 
了 先前 的 孔 板 模型 ， 使 得 压力 测试 可 用 于 计算 发 动 
机 内 的 封 严 效率 , 并 使 用 CFD 的 方法 与 实验 数据 验 
证 . Lalwani 等 03] 进行 了 定常 计算 ， 比 较 了 不 同 冷 
气量 下 轮 缘 密 封 盘 腔 内 的 流动 及 封 严 特性 。 

由 于 现今 预 旋 盘 腔 的 研究 中 通常 将 轮 缘 密封 结 
构 简 化 , 分 开 研究 预 旋 盘 腔 与 轮 缘 密封 . 本 文 将 进 气 
预 旋 与 轮 缘 密封 结合 在 一 起 ， 研 究 进 气 预 旋 与 轮 缘 
密封 的 相互 影响 . 本 文 以 文献 [14 中 模型 为 基础 , 采 
用 轴 向 进 气 ， 对 带 径 向 轮 缘 密封 的 预 旋 盘 腔 系统 进 
行 CFD 数值 研究 , 重点 研究 在 旋转 雷诺 数 0.8 x 108 
下 ,不 同 无 量 纲 冷却 气体 流量 和 进 气 预 旋 比 的 范围 
内 ， 熏 腔 内 的 流动 换 热 特性 ， 深 入 分 析 进 气 预 旋 对 
封 严 效率 的 影响 和 考虑 轮 缘 密封 后 盘 腔 内 部 流动 传 
热 特性 的 变化 ， 


1 计算 模型 和 数值 方法 


图 1(a) 给 出 了 计算 模型 静止 域 、 旋 转 域 、 动 静 
交界 面 及 盘 腔 中 分 面 的 位 置 示意 。 根据 文献 [15] 中 
的 研究 ， 将 动静 交界 面 设 定 在 盘 腔 左 侧 ， 静 叶 流 道 
计算 域 为 静止 域 ,动静 交界 面 右 侧 的 动 叶 流 道 和 盘 
腔 计算 域 为 旋转 域 . 盘 腔 左 侧 壁面 为 静止 壁面 , 右 侧 
壁面 随 旋转 域 旋转 . 同时 图 1(b) 给 出 了 计算 模型 的 
结构 示意 图 ， 其 主要 参数 如 表 1 所 示 . 冷却 空气 从 
盘 腔 左下 角 的 进 气 口 通 入 盘 腔 , 并 带 有 一 定 预 旋 度 ， 
对 转盘 壁面 进行 冷却 ， 最 后 经 轮 缘 处 的 径 向 密封 进 
入 主流 ， 起 到 封 严 作用 。 同时 主流 燃气 经 轮 缘 间 阶 
入 侵 盘 腔 。 盘 腔 内 的 流动 换 热 同时 受 冷却 气体 和 入 
侵 燃 气 的 影响 ， 

采用 商用 软件 进行 模型 构建 网 格 划分 和 数值 计 
算 . 图 2 为 运用 软件 NUMECA 生成 的 多 块 结构 化 计 
算 网 格 , 主流 通道 采用 H-O-H 拓扑 结构 化 网 格 . 主 
流通 道中 静 叶 周 向 节点 为 49 个 , 轴 向 节点 为 49 个 ; 
动 叶 周 向 节点 为 67 个 , 轴 向 节点 为 33 个 ; 主流 沿 叶 
高 方向 布置 了 37 个 节点 ; 为 保证 求解 精度 ， 生 成 的 
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主流 通道 与 轮 缘 密封 间隙 处 网 格 结 点 完全 匹配 ,在 
间隙 处 布置 了 37 个 网 格 节点 ; 网 格 最 小 正 交角 度 为 
32°, 质量 良好 。 节点 总 数 为 162 万 , 满足 网 格 无 关 
性 要 求 。 主流 计算 域 壁面 第 一 层 厚度 设 为 0.01 mm， 
盘 腔 计算 域 壁面 第 一 层 厚度 为 0.001 mm， 以 满足 采 
用 软件 ANSYS-CFX 数值 求解 三 维 RANS 方程 进行 
数值 计算 研究 时 ，SST 汕 流 模型 对 壁面 网 格 yt 的 
要 求 . 


图 1 主流 盘 腔 结构 与 计算 域 模型 示意 图 


Fig. 1 The mainstream and cavity configuration and 


computational domain 


表 1 径 向 轮 缘 密 封 预 旋 盘 腔 的 主要 结构 参数 (mm) 
Table 1 Dimensions of the pre-swirl rotor-stator 


cavity with rim radial seal(mm) 


结构 参数 数值 
内 半径 a 120 
外 半径 190 
间距 c 5 
动静 叶 间距 d 12 
冷气 进 气 位 置 rp 121 
主流 道 径 向 高 度 h 10 
盘 腔 间距 s 20 
冷却 气体 进口 宽度 sin 2 
盘 腔 轴 向 间隙 sax 5.3 
轮 缘 密 封 轴 向 间隙 sc,ax 2 
轮 缘 密 封 径 向 间隙 sc rad 2.4 


流体 介质 采用 理想 气体 , 转速 3000 r/min; 主流 
进 气 给 定 进 气量 0.51 kg/s 和 总 温 320 K; 冷气 进口 
给 定 进 气 速度 矢量 和 静止 参考 系 下 总 温 300 K; 主 
流出 口 以 大 气压 强 为 静止 参考 系 平均 静 压 ; 固 壁 面 
均 设 为 无 滑 移 壁 面 ; 盘 腔 静止 壁面 绝热 ;转动 壁面 
采用 给 定 温度 293 K 与 绝热 壁面 两 种 边界 条 件 ， 以 
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便 Nw 的 分 析 研 究 ; 计算 中 动静 交界 面 以 冻结 转子 
法 (Frozen Rotor) 连接 ; 两 侧 轮 盘 周期 面 设 为 旋转 
周期 关联 面 . 当 连 续 方 程 、 能 量 方程 、 满 流 方程 和 动 
量 方程 的 均 方 根 残 差 均 小 于 10-5， 且 观测 点 附加 变 
量 值 保持 不 变 时 , 认为 计算 收敛 . 


图 2 径 向 轮 缘 密封 盘 腔 模 型 的 计算 网 格 


Fig. 2 Computational mesh of the pre-swirl rotor-stator 


cavity with rim radial seal 


实验 中 测量 盘 腔 轮 缘 密封 的 封 严 效率 采用 浓度 
法 ， 即 在 主流 中 加 入 一 定 浓 度 的 示 踪 气体 ， 通 过 盘 
腔 内 固定 测 点 示 踪 气体 浓度 表征 轮 缘 密封 的 封 严 效 
率 . 封 严 效 率 se 定义 为 : 

sc 一 (cs 一 caj/(co 一 ca) (1) 
式 中 ,cs 为 主流 进 气 时 示 踪 气体 的 浓度 ，c 为 参考 
点 的 示 踪 气体 浓度 ，oo 为 冷气 流 进口 处 示 踪 气体 的 
浓度 . 

本 文采 用 添加 附加 变量 的 方法 模拟 示 踪 气体 在 
主流 通道 和 盘 腔 内 部 的 输 运 扩散 ,冷气 进口 示 踪 变 
量 CO, 设 为 1, 主流 进口 设 为 0, 并 设置 动能 扩散 系 
数 Dy 辅助 计算 (300 K 标准 大 气压 下 , COs 在 空气 
中 Dy = 1.6 x 10 m/s). 由 于 冷气 进口 示 踪 变 量 
CO 为 1, 主流 进口 为 0, 则 co=1，ca = 0， 封 严 效 
率 se 等 于 附加 变量 COs 的 值 . 


2 结果 与 讨论 
2.1 盘 腔 流动 特性 

盘 腔 内 预 施 比分 布 直接 反映 进 气 预 施 对 流动 的 
影响 . 定义 进 气 预 旋 比 Bu， 


bp = Veo,p/Tp! (2) 


式 中 , Vp 为 进 气 位 置 中 分 面 处 气流 切 向 速度 ; 0 为 
转盘 转速 。 
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定义 水 流 系数 XT: 
MAT = CuwRez0s (3) 


图 3 给 出 了 不 同 冷却 气体 进 气 量 , 不 同 进 气 预 
旋 比 时 盘 腔 中 分 面 >/s=0.5 处 预 旋 比 沿 径 向 的 分 布 . 
比较 图 3(a)~(c)， 进 气 预 旋 比 越 大 盘 腔 内 预 旋 比 越 
大 ; 同时 , 随 着 冷气 进 气 量 逐 渐 增 大 , 相应 的 进 气 预 
旋 比 变化 对 盘 腔 内 流动 影响 加 大 . 

当 Cw = 1500 时 , 和 Tr = 0.028, 进入 盘 腔 的 冷气 
流 受 转 静 盘 壁 面 边界 层 影 响 极 大 ， 流 动 处 于 秋 性 主 
导 ; 随 着 径 向 位 置 的 增 大 , 进 气 流 与 盘 腔 原 有 气流 逐 
渐 混 合 , 不 同 不下 的 8 分 布 受 盘 腔 内 气流 影响 相应 
增 大 减 小 , 直至 7/b=0.9 时 趋 于 一 致 . 当 Cw = 3500 
时 ，AT = 0.069， 流 动 逐步 从 秋 性 主导 向 惯性 主导 
过 度 ，9 对 盘 腔 内 预 旋 比 分 布 影响 增 大 ,代表 不 同 
Bop 的 四 条 曲线 更 加 分 散 ， 中 分 面 旋转 核心 区 受 壁 面 
影响 减 小 . 当 Cw = 6000 时, 和 rT = 0.113, 流动 为 惯性 


0.65 0.70 0.75 0.80 0.85 
r/b 


(a) Cw=1500 


0.90 0.95 


0.0 


0.65 0.70 0.75 


0.80 0.85 
r/b 
(b) Cw=3500 


0.90 0.95 


ChinaXiv 合 作 期 
1649 


平 | 


| 


201711.02573V1 


chinaXiv 


1650 工程 热 物 理学 报 


0.80 0.85 
r/b 


(c) Cw=6000 
图 3 预 旋 比 沿 径 向 分 布 


Fig. 3 Radial variation of pre-swirl ratio 
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主导 , 漠 动 剧烈 , 7/b = 0.64 处 由 于 进 气流 的 影响 可 
以 明显 看 到 6 波动 ; 由 于 壁面 边界 层 影 响 减 小 , 至 
7 三 0.95 处 0 才 急 剧 增 大 . 
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(b) Cw=3500 
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(c) Cw=6000 
图 4 r = 182.6 mm 处 预 旋 比 轴 向 变化 


Fig. 4 Axial variation of pre-swirl ratio at 7 = 182.6 mm 


图 4 为 不 同 冷却 气体 进 气 量 , 不 同 进 气 预 旋 比 
时 盘 腔 气流 进入 密封 时 (7 = 182.6 mm) 预 旋 比 沿 轴 
向 的 分 布 . 当 Cw = 1500, Ar = 0.028 时 ,如 图 4(a) 
所 示 , 流动 处 于 茜 性 主导 , 沸 动 较 小 , 进 气 与 盘 腔 内 
气体 混合 均 句 ,气流 流 至 密封 时 已 不 受 进 气 预 旋 的 
影响 , 此 时 转 静 壁面 边界 层 主导 盘 腔 内 流动 状态 . 当 
进 气流 逐渐 增 大 , 消 动 增强 , 如 图 4(b)、(c) 所 示 , 流 
动 受 边 界 层 影 响 逐 渐 减 弱 ， 进 入 密封 时 气流 仍 受 进 
气 预 旋 影 响 , 研究 密封 处 流动 时 思 影响 不 可 忽略 . 
同时 随 着 近 壁 面 边界 层 影 响 减 小 ， 进 入 密封 气流 的 
6 整体 减 小 . 

图 5 为 不 同 冷却 气体 进 气 量 ， 不 同 预 旋 比 时 
预 旋 比 在 盘 腔 内 子午 面 的 分 布 云图 。 如 图 所 示 ， 当 
Cw = 1500 时 ， 进 气 预 旋 比 变化 仅 影响 到 盘 腔 中 上 
部 , 盘 腔 顶部 9 分 布 相似 . Cw = 3500 时 , 进 气流 系 
数 对 盘 腔 整 体 流动 影响 加 强 . Cw = 6000 时 , 进 气 预 
旋 明 显影 响 密 封 内 流动 状态 。 


Cw=6000 


Cw=1500 


Cw=3500 


图 5 预 旋 比 在 盘 腔 内 子午 面 的 分 布 云图 


Fig. 5 Distribution of pre-swirl ratio in the meridional plane 


定义 无 量 纲 静 压 系数 &: 
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§€=2(p— pp)/(p2°b) (4) 


式 中 , pp 为 中 分 面 进 气 处 静 压 . 

图 6 为 不 同 冷却 气体 进 气 量 ,， 不 同 进 气 预 旋 比 
时 静 压 系 数 沿 径 向 的 变化 . 如 图 6 所 示 , 由 于 旋转 泵 
效应 , 静 压 随 着 半径 的 增 大 存在 明显 的 压力 梯度 , & 
沿 径 向 线性 增 大 , 且 其 斜率 主要 受 进 气 Pp 影响 . 图 
6(a) 中 进 气流 量 较 小 ， 盘 腔 内 流动 受 转 静 盘 边 界 层 
影响 大 ， 旋 转 泵 效应 愈加 明显 ， 曲 线 斜率 较 大 . 图 
6(b), (c) 中 , 流动 处 于 惯性 主导 , 受 流量 影响 较 小 ， 
斜率 由 进 气 Bb 主导 . 此 外 , 当 Bb = 1 时 , 《 分布 几 
乎 不 受 流量 影响 。 

轮 缘 密 封 处 燃气 入 侵 现象 出 现在 主流 压力 高 于 
盘 腔 压力 处 ， 进 气 预 旋 增 大 盘 腔 静 压 ， 由 此 推测 进 
气 预 旋 有 利于 轮 缘 密 封 封 严 效率 的 提高 ， 但 在 冷气 
流量 较 小 时 ， 受 转 静 盘 边界 层 影响 ， 进 气 预 旋 对 盘 
腔 内 流动 影响 较 小 ， 静 压 变 化 不 大 ， 因 此 封 严 效率 
提高 有 限 . 


明寺 可 五 
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(c) Cw=6000 
图 6 静 压 系数 沿 径 向 变化 


Fig. 6 Radial variation of static pressure coefficient 
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定义 无 量 纲 总 压 系数 60: 
é€0 = 2(po — po,p)/ (p220°) (5) 


式 中 po,p 为 中 分 面 进 气 处 总 压 . 
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(c) Cw=6000 
图 7 总 压 系数 沿 径 向 变化 


Fig. 7 Radial variation of total pressure coefficient 
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图 7 为 不 同 冷却 气体 进 气 量 ， 不 同 进 气 预 旋 比 
时 总 压 系数 沿 径 向 的 变化 . 与 上 变化 不 同 的 是 ， 盘 
腔 内 6 的 沿 径 向 基本 不 变 ; 由 总 压 定义 可 知 ,， 总 压 
与 静 压 ， 与 速度 平方 成 正比 ， 由 于 盘 腔 环境 中 切 向 
速度 较 大 , 因此 6 与 5 和 82 成 正比 , 符合 如 图 3、 
图 6 中 6 沿 径 向 减 小 , & 沿 径 向 增 大 的 趋势 . 对 比 图 
7(a)~(c),， 当 冷气 量 逐 渐 增 加 时 , 流动 受 近 壁面 边界 
层 影 响 减 小 , 单位 冷气 量 受 转盘 壁面 作用 功 减 小 , 6 
减 小 . 同时 , 随 着 进 气 预 旋 比 的 增 大 , 相对 转盘 气流 
速度 减 小 , 转盘 对 气流 的 作用 功 减 小 , So 减 小 , 级 效 
率 提高 . 此 外 , 在 进 气 位 置 附近 , 进 气流 与 盘 腔 内 气 
流 混 合 , 并 冲击 到 转盘 表面 , 一 部 分 动能 转化 为 热 ， 
因此 总 压 急剧 下 降 ， 尤其 是 当 fb, = 1.0 时 , 进 气 初 
速度 较 大 , 混合 冲击 均 较 剧烈 ,6o 明显 较 小 ; 而 高 位 
近 上 壁 处 ， 受 壁面 影响 ， 切 向 速度 迅速 增加 ， 和 co 上 
升 : 
2.2 盘 腔 换 热 特性 

表征 换 热能 力 大 小 的 努 寨 尔 数 Nu 定义 如 下 : 


Nu = hr/k = qwrT/ [k(Tw 一 Twad)] (6) 


式 中 , h 为 表面 传 热 系数 ; qw 为 壁面 热流 密度 ;上 为 
导热 系数 ; Tvw 为 壁面 温度 ; Tw,ad 为 壁面 绝热 温度 . 

在 过 去 的 研究 中 ， 转 静 盘 系统 壁面 绝热 温度 
Twad 常 依据 文献 [16]， 由 下 式 根据 盘 腔 流动 计算 得 
出 : 


0 人 2 vp,00 \? 
re (1- 区) 
I Br 


式 中 ,Th 。 为 进 气 总 温 ， vw, 为 旋转 核心 处 切 向 速 
度 ; RR 为 温度 恢复 系数 ( 取 Pr1/3). 


由 于 轮 缘 处 存在 燃气 入 侵 , 为 研究 入 侵 燃 气 对 
盘 腔 内 部 的 换 热 的 影响 ,比较 式 (6) 与 式 (7) 计算 得 
出 的 努 塞 尔 数 结果 . 图 8 为 不 同 冷却 气体 进 气量 , 不 
同 预 旋 比 时 Nw 沿 径 向 分 布 . 如 图 所 示 , 当 Ow 较 小 ， 
燃气 入 侵 较 严重 时 ,Nu 明显 大 于 Nuz, 而 当 Cw 较 
大 , 封 严 较 好 时 ， 两 种 方式 计算 所 得 Nu 相差 不 大 ; 
因此 可 知 当 存在 燃气 入 侵 时 ， 由 于 燃气 较 高 温度 的 
影响 , 盘 腔 内 气流 温度 受热 上 升 , 实际 Tw,aa 大 于 根 
据 旋转 核心 区 流动 状态 计算 所 得 , Nuz 小 于 Nu; 而 
当 封 严 较 好 时 ，Nw 与 Nu 相近， 盘 腔 内 换 热 计算 
可 忽视 主流 入 侵 影 响 . 此 外 , 进 气 b, 的 增 大 会 导致 
转盘 进 气 高 度 处 Nu 局 部 增 大 , 这 种 Nw 的 局 部 变 
化 两 种 计算 方法 均 可 得 到 , 且 随 着 冷气 量 的 增加 , 滑 
动 的 增强 , 变化 愈加 明显 . 


一- 量 : 一 :一 p=0, Nul 
= b=0, Nu, 
ee f=0.5, Nu 


-一 ee- p=0.5, Nu 
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9 0.65 0.70 0.75 0.80 0.85 0.90 0.95 
r/b 
(b) Cw=3500 
(Nu 一 一 根据 式 (7) 所 得 ，Nus 一 一 根据 定义 式 (6) 所 得 ) 


图 8 转盘 表面 Nw 沿 径 向 分 布 


Fig. 8 Radial variation of the nusselt number on rotor 
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图 9 中 无 量 纲 温度 系 数 Ow 定义 如 下 : 
Ow = (Tw,aq To,p)/ (To,in 一 70p) (8) 


式 中 To,in 为 主流 进 气 总 温 . 

图 9 为 不 同 冷却 气体 进 气量 ,不 同 预 旋 比 时 转 
盘 壁 面 温度 变化 云图 . 在 轮 缘 处 , 入 侵 燃 气 对 转盘 的 
加 热 主 导 壁 面 温度 的 变化 . 如 图 9 所 示 ， 入 侵 燃 气 
对 转盘 密封 处 壁面 的 绝热 温度 影响 较 大 ， 且 存在 明 
显 的 周 向 不 均匀 性 . 当 Cw = 1500 时 , 封 严 较 差 , 入 
侵 量 大 ， 将 影响 整个 盘 腔 的 温度 分 布 ; 随 着 冷气 量 
的 逐步 增加 ， 封 严 效 率 的 提高 ， 入 侵 气 流 对 密封 处 
转盘 壁面 及 盘 腔 内 的 温度 影响 也 逐渐 减 小 。 对比 不 
同 进 气 预 旋 比 时 的 温度 云图 变化 , 当 Bb = 1.0 时 入 
侵 对 密封 壁面 温度 影响 减弱 , 但 局 部 温度 梯度 增加 ; 
当 进 气 预 旋 增 大 时 ， 入 侵 气 流 对 密封 壁面 的 影响 控 
制 在 了 更 小 的 范围 内 . 


| 
Cw=1500 Cw=3500 
图 9 转盘 壁面 温度 云图 


Fig. 9 Distribution of adiabatic temperature on rotor surface 


2.3 轮 缘 密 封 封 严 特 性 

封 严 效率 sc 由 式 (1) 定义 ， 取 静 盘 表面 "一 
182 mm 处 se 决定 径 向 密封 的 封 严 效率 . 图 10 是 
3 种 进 气 预 旋 比 下 轮 缘 密封 的 封 严 效率 随 冷 气量 变 
化 折线 图 ， 如 图 所 示 ，B。 的 增加 提高 径 向 密封 封 
严 效 率 , 但 在 冷气 量 较 小 时 影响 有 限 ; 冷气 量 达到 
Ow = 6000 时 , ec 达到 1. 

主流 压力 系数 所 定义 为 : 


ka = 2(pa — Ba)/ (p22b?) (9) 


式 中 , ps 为 轮 载 处 主流 静 叶 尾 缘 与 轮 缘 间 隙 中 间 位 
置 处 的 静 压 ; 5。 为 pa 的 平均 值 . 

图 11 为 主流 通道 处 周 向 压力 分 布 ， 如 图 所 示 ， 
主流 存在 周 向 的 压力 波动 ， 波 动 导致 主流 气体 对 盘 
腔 的 外 部 诱导 入 侵 (ET), 这 种 入 侵占 据 主 流入 侵 的 
主要 部 分 台 ; 主流 压力 受 冷 气量 与 冷气 进 气 预 旋 比 
的 影响 不 大 . 
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图 10 不 同 冷气 量 下 轮 缘 密 封 的 封 严 效 率 


Fig. 10 Effect of sealing flow rate on sealing efficiency 
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图 11 主流 通道 处 周 向 压力 分 布 
Fig. 11 Circumferential distribution of static pressure 
coefficient in annulus 


图 12 为 密封 间隙 处 径 向 速度 云图 . 如 图 12 所 
示 ， 即 使 冷气 量 达到 封 严 所 需 ， 主 流 气 流 由 于 高 压 
作用 仍 会 对 密封 间 队 有 一 定 的 入 侵 ; 主流 气流 在 周 
向 和 峰值 出 入 侵 盘 腔 ， 与 封 严 气流 混合 ,一 部 分 流 
入 盘 腔 内 部 造成 更 深度 的 主流 入 侵 ， 大 部 分 随 封 严 
气流 自 & 谷 值 位 置 附近 流出 间 际 . 

图 13 为 径 向 密封 不 同 高 度 在 7 = 182.6 mm， 
188.8 mm 和 194 mm 处 的 封 严 效率 云图 ， 如 图 
13 所 示 ， 随 着 冷气 量 的 增加 ， 主 流入 侵 的 深度 减 
小 . Cw = 1500 时 ， 沸 动 较 小 ， 气 流 混合 均匀 ， 密 
封 处 周 向 不 均匀 性 较 低 ,同时 由 于 混合 充分 , 在 & 
较 小 的 位 置 入 侵 气 流 无 法 排 尽 ， 导 致 盘 腔 的 封 严 效 
率 较 低 . 当 Cw = 6000 时 , 达到 最 小 封 严 气量 , 主流 
的 入 侵 对 盘 腔 的 影响 基本 控制 在 密封 的 轴 向 间隙 密 
封 段 . 
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图 12 轴 向 密封 间 际 处 径 向 流速 云图 


Fig. 12 Distribution of radial velocity in axial clearance 
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图 13 径 向 密封 处 封 严 效率 多 层 云 图 


Fig. 13 Distribution of sealing efficiency at different radius 


3 结 论 


本 文采 用 附加 变量 法 ,数值 求解 了 三 维 RANS 
方程 组 和 SST 清流 模型 , 分 析 了 三 种 不 同 冷 气量 下 ， 
四 种 冷气 进 气 预 旋 比 时 ， 带 径 向 轮 缘 密封 和 轴 向 预 
旋 进 气 的 转 静 盘 腔 的 流动 换 热 特 性 与 轮 缘 处 密封 的 
封 严 特性 , 得 出 如 下 结论 

1) 盘 腔 内 中 分 面 处 预 旋 比 随 半径 的 增 大 而 减 小 ， 
其 总 体 随 进 气 预 旋 比 的 增 大 而 增 大 ; 中 分 面 处 静 压 
系数 沿 径 向 线性 增 大 ， 和 斜率 大 小 与 进 气 预 旋 比 成 正 
比 ; 中 分 面 处 总 压 系数 与 预 旋 比 平方 及 静 压 系数 成 
正比 , 沿 径 向 基本 不 变 , 略 有 增 大 . 

2) 当 入 侵 较 大 时 转盘 壁面 绝热 温度 升 高 ， 壁 面 
绝热 温度 无 法 通过 理论 公式 由 盘 腔 内 流动 状态 得 出 ， 

理论 式 得 出 的 Nu 较 实 际 略 大 . 此 外 主流 入 侵 对 
密封 处 壁面 温度 影响 较 大 ， 使 壁面 绝热 温度 出 现 明 
显 的 周 向 温度 梯度 ， 温 度 梯度 随 进 气 预 施 比 增 大 而 
增 大 . 

3) 封 严 效率 随 冷气 量 增 大 而 增 大 , 进 气 预 旋 提 
高 封 严 效率 . 随 冷 气量 增加 , 入 侵 气 流 与 盘 腔 内 的 气 
流 掺 混 控 制 在 密封 轴 向 间隙 与 小 部 分 径 向 间隙 中 . 
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